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Abstract 

A method for finding relations for the roots of polynomials is presented. 
Our approach allows us to get a number of relations for the zeros of the 
classical polynomials and for the roots of special polynomials associated 
with rational solutions of the Painleve equations. We apply the method to 
obtain the relations for the zeros of several polynomials. They are: the La- 
guerre polynomials, the Yablonskii - Vorob'ev polynomials, the Umemura 
polynomials, the Ohyama polynomials, the generalized Okamoto polyno- 
mials, and the generalized Hermite polynomials. All the relations found 
can be considered as analogues of generalized Stieltjes relations. 

Keywords: Special polynomials, Rational solutions, Stieltjes relations, the 
Painleve equation, power expansion 

PACS: 02.30.Hq - Ordinary differential equations 

1 Introduction 

In 1885 Stieltjes [1, 2] suggested the following relations for the roots of the 
Hermite polynomials [3,4] 



^2(Zj-z k ) 1 -z j = 0, (j = l,...,n) (l.l) 

These relations can be considered as the equilibrium configurations of n 
particles on the line interecting pairwise with a repulse logarithmic potential in 
the harmonic field [4-6]. 

In this paper we introduce a new approach for finding relations like (|1.1|) 
satisfied by the roots of polynomials, which can be generated by a differential 
equation. We demonstrate our approach to get relations for the roots of the 
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Laguerre polynomials and for the roots of special polynomials associated with 
rational solutions of the second, the third, and the fourth Painleve equations. 

The six Painleve equations (Pi — Pq) were picked out by Painleve and his 
student Gambier from a certain class of second - order differential equations as 
equations whose general solutions on the one hand, do not have movable critical 
points and on the other hand, cannot be expressed through known elementary 
or special functions. In other words they define new transcendental functions. 
However, the equations P2 — Pe possess hierarchies of rational and algebraic 
solutions at certain values of the parameters. It turned out that such solutions 
could be described using families of special polynomials. The history of this 
question is as follows. 

A. I. Yablonskii and A. P. Vorob'ev were the first who expressed the rational 
solutions of P 2 via the logarithmic derivative of the polynomials, which now 
go under the name of the Yablonskii - Vorob'ev polynomials [7,8]. Later K. 
Okamoto suggested special polynomials for certain rational solutions of P4 [9]. 
H. Umemura derived analogous polynomials for some rational and algebraic so- 
lutions of P3 and P5 [10]. All these polynomials possess a number of interesting 
properties. For example, they can be expressed in terms of Schur polynomi- 
als. Besides that the polynomials arise as the so - called tau-functions and 
satisfy recurrence relations of Toda type. Recently these polynomials have been 
intensively studied [11—13]. 

Apart from rational and algebraic solutions P2 — Pq possess one - param- 
eter families of solutions expressible in terms of the classical special functions. 
In particular cases these special function solutions are reduced to classical or- 
thogonal polynomials and, consequently give rational solutions of the Painleve 
equation in question. For P4 these polynomials arc Hermitc polynomials, for P3 
and P5 Laguerre polynomials and for Pg Jacobi polynomials. 

Not long ago P. A. Clarkson and E. L. Mansfield suggested special polynomi- 
als for the equations of the P2 hierarchy [14]. Also they studied the location of 
their roots in the complex plane and showed that the roots have a very regular 
structure. 

The textbook [15] characterizes the Painleve equations as "the most impor- 
tant nonlinear ordinary differential equations" and states that " many specialists 
believe that during the twenty-first century the Painleve equations will become 
new members of the community of special functions" [16]. 

In this paper we study the relations for the zeros of the following equations 



yijzzzz 



- ^VzVzzz + 3y zz 2 - 



z (yy Z z - Vz 2 ) -yy z = 



(1.2) 



z 2 (y Vzzzz - 4 y z Vzzz + 3 y 2 zz ) + 2 z (y y zzz - y z y zz ) 
-4z(z + [i)(yyzz - vl) - 2yy zz + 4nyy z - 2n(n+ 1) y 2 = 0, 



(1.3) 



z 3 {y y ZZ zz -4y 2 y Z zz + iy 2 Z z) ~ 6z 2 {yy ZZ z - Vz yzz) 
-12 z{z i - 3n-l)(yy zz -y 2 z )- 9 z (y y zz + y 2 z ) 
+3(12 z 4 - 16 n z 2 + 12 n + 7)y y z 
-24nz((n + 3)z 2 -3n-l)y 2 = 0, 



(1.4) 
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.'///-.:-. . - ^VzVzzz + 5y 2 zz + ^z 2 (yy zz - y 2 z ) + izyy z 

g 

— (m 2 + ft + ran — m — n) y 2 = 0, 

o 



(1.5) 



y Vzzzz ~^y z Vzzz + 3y 2 zz ~ ^(z 2 + 2n - 2m) {yy zz ~y 2 z ) . * 

+Azyy z — 8mny =0. 

It is known [11] that particular solutions of equation i|1.2|) are the Yablonskii 
- Vorob'ev polynomials. Rational solutions of the second Painleve equation 
can be expressed as the logarithmic derivative of these polynomials. In this 
work using our approach we obtain some new relations for the zeros of these 
polynomials. These relations explain a very regular location of the zeros in the 
complex plain. 

The equations l|1.3|) and (|1.4|) in their turn possess polynomial solutions, 
i.e. the Umemura polynomials and the Ohyama polynomials [12], accordingly. 
Rational solutions of the third Painleve equation can be obtained with the help 
of these polynomials. The general form of these polynomials is not known at 
present but using the method of our recent works [17, 18] one can solve this 
problem and can get the general form for these polynomials. However here we 
present new additional correlations for the zeros of these polynomials. 

The equations (|1.5(l and (|1.6(l also admit polynomial solutions, which are 
the generalized Okamoto polynomials and the generalized Hermite polynomials, 
respectively [13]. Rational solutions of the fourth Painleve equation can be 
expressed via logarithmic derivative of these polynomials. The general form of 
these polynomials is not known as well. In this work new correlations for the 
zeros of these polynomials are obtained. 

This paper is organized as follows. In section 2 we present the method 
of finding the relations for the roots of the polynomials. In this section we 
also apply our approach on the example to have the relations for the roots 
of the Laguerre polynomials. Sections 2, 3, 4 and 5 are devoted to finding 
the relations for the zeros of the Yablonskii - Vorob'ev polynomials, of the 
Umemura polynomials, of the Ohyama polynomials, of the generalized Okamoto 
polynomials and of the generalized Hermite polynomials. All these polynomials 
are associated with rational solutions of the second, the third and the fourth 
Painleve equations. 



2 Method applied 

Let us assume that we have an n-linear (linear, bilinear etc) equation of order 



E 1 (y,y z ,y zz ,...,z) = 0. (2.1) 

In addition suppose that this equation admits polynomial solution y(z) = q{z) 
of degree p. The transformation 

w(z) = V — (2.2) 

y 
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maps solutions of the equation l|2.1f> into solutions of the nonlinear equation 

E2(w,w z ,w zz ,..., z) = 0. (2.3) 

Note that the degree of 1)2. 3|l is M — 1. Under such transformation the polyno- 
mial q{z) becomes the rational function 

^) = E-^, (2-4) 

where Zj is a root of q(z) of degree ej and 
p 

p = ]Te„ (2.5) 

i=i 

i.e. p is the number of distinct roots. The rational function (|2.4|l can be ex- 
panded in a neighborhood of a pole Zj and this gives 

oo / 



< I z — Zj | < min | z^ — Zj | 

At the same time we can find expansions of solution w{z) near movable 
singular point zo using the Painleve test [19-22] or methods of power geometry 
[23-27]. The result of calculations for the rational function w(z) takes the form 



(z)= 5> m (z-z r-\ (2.7) 

m=0 

Assuming zq = Zj, cq — ej and comparing this series with l|2.6|) yields 



E (zfe _ e ")m + 1=- C - m >°- (2.8) 

We are going to use correlations 1)2. 8JI for finding some relations of the zeros 
of the above mentioned polynomials, which are associated with the rational 
solutions of the second, the third, and the fourth Painleve equations. 

Let us demonstrate our approach by the example of the Laguerre polynomials 
Ln*\z). They satisfy the following second-order linear differential equation 

zy zz + {a + l-z)y z + ny = 0. (2.9) 

Substituting y z = w(z) y{z), y zz = (w z + w 2 )y into (|2.9(l and dividing the result 
by y(z) we obtain the first-order nonlinear differential of the form 

zw z + zw 2 + (a + 1 - z)w + n = 0. (2.10) 
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Taking the Painlcvc methods into account we have the power expansions for 
rational solution of equation l|2.10[l near singular point z = zq (zq ^ 0) in the 
form 



1 1 a + 1 - Z(] 1 
w(z) = 1 =■ 6 a + 5 — 2 zq — 2 a zq- 

Z — Zq 2 Zq 12 Zq 



+zq 2 — 4 uzq + a 2 ] (z — Zq) - [ 4 a + a 2 + 3 — a zq — 2 uzq- 

8zo 3 



-zq } (z - z f - 4 [ 76 a z - 360 a - 8 uzq a 2 + 2 z 2 + ^ ^ 



-152 nz - 110 a 2 -4a 2 z + 8az Q 2 + 16nz 2 - 4z 3 - 251- 
+z 4 + a 4 - 8 z 3 n + 6 a 2 z 2 - 4 a z 3 - 4 a 3 z + 



+ 16 n 2 zo 2 + 16 a ZQ 2 n + 76 zq ] (z — zq) 3 + ... 

Comparing the power expansion (|2.11l) with 12.6(1 we have the following 
relations for the zeros of the Laguerre polynomials in the case Zj 



2^(z k - Zj ) 2 Zj 2 



~=0, (2.12) 



efe a 2 + 6 a + 5 1 + a + n 1 



^(z fe -z^ + 12 z 2 6z 3 +12 °' (2 ' 13) 



8z 3 



1 + a + 2 n 
8^f 



= 



(2.14) 



In this expressions n makes sense of p and e k = 1, Zf. ^ 0. It is known that 
z = is the root of L^^z) if and only if a = — I (1 < I < n). Then in this case 
the asymptotic series in a neighborhood of z = is 



w(z) 
(5n 2 l + 



_ I 

z 

11 n 2 



l + l 
- lln 



I (n~l)(n+l)z 
(1 + 2) (l + l) 2 
5 nl 2 - 6 nl + 2 - 



(2n+l-l)z 
_ 1+ (1 + 3) (l + l) 

6l + l 3 -l 2 )z 2 



(l + S)(l + 4)(l + 2)(l + iy 



+ o(z 2 ) 



(2.15) 



There are no arbitrary constants in this expansion as the critical number k 
is k = — I — 1 < — 1. Without loss of generality let us suppose that z = is the 
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first in the set of roots. Comparing the series 1)2. 15[) with l|2.6|) we see that the 
order of the root z — is I and 

h , 

El n — I , 
^—l- * (2 - 16) 



EA- ( "" 0( " + 1 No, (2-17) 
f^ 2 zl (l + 2)(l + lf 



^ _L (n-l){n+l){2n+l-l) 

z 



, — ^- - (2.18) 

(Z + 2)(/ + 3)(/ + l) 3 



We would like to note that in such a way one can find analogous relations 
for some other polynomials. 

3 Relations for the zeros of the Yablonskii - Vorob'ev 
polynomials 

The transformation H2.2(l maps solutions of the equation (|1.2|) into solutions of 
the third-order nonlinear differential equation 

w zzz + 6w z 2 — zw z — w = 0. (3-1) 

The Laurent expansion of its solution w{z) in a neighborhood of a movable pole 
zq is the following 

2 



w{z) = Jz^ +ai+ l2 iz - Zo) -{720 + 3d^ Z - Za) 



"144 {Z ~ Z0) + a ^ z ~^ " I 8640 + 360 ) + 



(3.2) 



The Yablonskii - Vorob'ev polynomials Q n (z) are monic polynomials of degree 
n(n + l)/2 with simple roots. Then from power expansion 1)3. 2|) we obtain the 
following relations for the zeros of Q n (z) 

n(n+l)/2 
k=£] 



n(n+l)/2 



n(n+l)/2 



E 



= 0, 



(3.5) 
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n(n+l)/2 n(n+l)/2 



E 777^y7 + ?n E T^TT-Sn- ' ^ 



n(n+l)/2 

E 

k^3 



E 7 — \E -Th=°> (37) 



n(n+l)/2 

n(n+l)/2 n(n+l)/2 2 

V 1 + _L V 1 - -^L_ =n f3 91 

^ (z k - z 3 y 360 £j (z k - zj) 8640 1 ' ; 

In this expressions the coefficients a\ and ag depend on n and Zj. Further 
the coefficient can be calculated if we note that the equation 13.1JI possesses 
the first integral 

W lz - \ W z* + 4w z _ ZW l - 2wUJ z + \ ZW + C = 0. (3.10) 

The latter equation has rational solutions in the form w(z; n) — d z \n(Q n (z)) if 
C = —n(n + l)/4. The Laurent expansion of w(z;n) = d z \n(Q„(z)) coincide 
with 13. 2[) accurate to the value of a§, which is no longer arbitrary 

1 3 1 1 1 1 2 fo-,^ 

dp, = Zi Zj ai H n H n . ( o.ll ) 

6 30240 3 280 420 1 560 560 y ' 

We do not present here the Yablonskii - Vorob'ev polynomials. They were 
obtained before [7,8, 11, 17]. Let us give here only the polynomial Q-r(z). That 



Q 7 (z) = [z 27 + 504 z 24 + 75 6 00 z 21 + 5174400 z 18 + 62092800 z 15 + 

+13039488000 z 12 - 828731904000 z 9 - (3.12) 
-49723914240000 z 6 - 3093932441600000] z. 

The configuration of the zeros for the Yablonskii - Vorob'ev polynomial Qt{z) 
is illustrated in the Figure 1. 



4 Relations for the zeros of the Umemura poly- 
nomials 

Using the formula y z = w(z) y{z) we can transform equation (|1.3|) to the fol- 
lowing third-order nonlinear equation 

z 2 w zzz + 6 z 2 w z 2 + 2 z w zz + 4 zww z — 4 z 2 w z — 4 z w z /i — 2 w z — 

(4.1) 

-2w 2 +4/iw-2n 2 -2n = 0. 
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Figure 1: Configuration of zeros for the Yablonskii - Vorob'ev polynomials Qy{z) 



With the help of the Painleve methods [19] we find the power expansion of its 
solutions w(z) near movable singularity zq 

w(z) = hai (z - z ) - 77 — q 4 «i z o~ 

Z — Zq 6 Zq 45Z() 

—2 zo 2 ai — 8 /i ai^o + 3 n 2 zo + 3 zq — 9 ai + 9 /i + 3 nzo + 2 /i zo 2 + 

(4.2) 

+^ 2 z + z 3 ] (z - z Q f + -^-j [fiz 2 - z 2 ai +2a 1 2 z - 4^iaiz + 

+[i 2 z + 3fi + n 2 z + z - 3ai + nz ] (z - z ) 4 + a 6 (z - z )° + ... 

Making use of the expansion near infinity for w(z) we easily find the degree 
p of the Umemura polynomials p = n(n + l)/2. Proceeding in the way stated 
in section 2 we obtain the following relations for the zeros of the Umemura 
polynomials in the case Zj ^ 



at = 0, 



(4.3) 



5(* 



^) 2 



3 Zj 3 



(4.4) 
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y 7 ^ yi = 0, (4.5) 



efc _ {(J, - Oi) _ 2 (/j - ai) _ J_ 
jL(z k - Zj )* 5z 3 45 Zj 45 

(4.6) 

(4ai 2 - 8^ai + /i 2 + 3 + 3n + 3n 2 ^ 



45^ 2 



0. 



Here e& = 1 for all Zfc 7^ 0. Substituting a\ from l|4.3|l into 1|4.4|) we obtain 

1 ^ 1 1 _|_ A* _|_ 1 Q (4 7) 

The Umemura polynomials Qs(z) takes the form [13] 



Q 5 (z) = z 15 - 135 z 14 + 8505 z 13 - 331380 z 12 + 8921745 z 11 - 

-175642425 z 10 + 2609560800 z 9 - 29764778400 z 8 + 
+262512608400 z 7 - 1788313917600 z 6 + 9322366561200 z 5 - (4.8) 
-36474213006000 z 4 + 103536532638000 z 3 - 200979491328000 z 2 + 
+238097325312000 z - 129459762432000 

The configuration of the zeros for the Umemura polynomial Q$(z) is demon- 
strated in the Figure 2. 

5 Relations for the zeros of the Ohyama poly- 
nomials 

Using y z = w(z) y we have from i|1.4|) the following third-order differential equa- 
tion 

z 3 w zzz -6z 2 w zz + 6z 3 w z 2 + (-12 z 2 w - 12 z 5 + 36 nz + 3 z) w z - 

-18 zw 2 + (21 + 36z 4 -48nz 2 + 36n)w + 24nz- (5.1) 
-24n 2 z 3 - 72 nz 3 + 72n 2 z = 0. 
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Figure 2: Configuration of zeros for the Umemura polynomial Qs(z) at /i = 16 

The power expansion for the solutions of equation 1)5. 1[) in a neighborhood of a 
pole z — zq (zq + 1 0) is 



w(z) = — h ai + — !-o (4 aiz - 7 - 12 n + 4 z 4 ) (z - z ) + 

Z ~ Zq 4 Zq 



+ 3 (21 - 16 aiz + 36 n + 20 z 4 - 16 ?iz 2 ) (z - z f 

8 z 



80z 



j [ 64 n 2 z 4 + 64 z^ ai z - 64 aiz 5 + 96 nz 4 + 120 z 4 + 144 1 



+128 naiz 3 + 16 z 8 + 301 - 320 nz 2 - 192 naiz - 192 n 2 z 2 - 



(5.2) 



-320 aiz + 600 n ] (z - z )' 5 + — [161 - 64 z^n + 64 n 2 z 4 

32z !5 



64z 2 ai 2 - 64aizo + 96nz 4 + 104 z 4 + 144 n z + 128naiz 3 



+16 z 8 - 208 nz 2 - 192naiz - 192n 2 z 2 - 208aiz + 



+360 n] (z - z ) 4 + a 6 (z - z ) 5 + 
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The degree of the Ohyama polynomial can be obtained with the help of the 
expansion of w(z) near infinity and equals n(n + 3)/2. From power expansion 
(|5.3I) we obtain the following relations for the roots of the Ohyama polynomials 
in the case Zj ^ 



^2 3 + oi = 0, (5.3) 



V- e fe 5 2n 2oi , (21 + 36n) 

g(^F + 2 z ^--T + ^^ = ' (5 - 5) 



e fc ^ 4a i (8n 2 + 12n+15) 8nai 



. ,. (z fe - zj) 4 5 5 J 10 5z 



(4a 2 -20n- 12n 2 ) 4ai(3n + 4) (144 n 2 + 301 + 600 n) 
5z] + bz] 80z/ 



(5.6) 



= 



In this expressions eu = 1 if Zfc ^ 0. Substituting ai from (|5.3(l into (|5.4|l we 
have 



i£j (*k - *i) 2 J 4z/ 

The Ohyama polynomial -Pz(z) can be written in the form [13] 



P 7 (z) = z 11 [z 24 - 56 z 22 + 1470 z 20 - 23 8 00 z 18 + 263375 z 16 - 
-2088240 z 14 + 12105940 z 12 - 51466800 z 10 + 158533375 z 8 - (5.8) 
-343343000 z 6 + 493643150 z 4 - 421821400 z 2 + 163788625] 

The configuration of the zeros for the Ohyama polynomial Pr{z) is presented 
in the Figure 3. 



6 Relations for the zeros of the generalized Okamoto 
polynomials 

Substituting the following expression y z — w(z)y into l|1.5|l sufficient amount of 
times we get 

w Z zz +6w z + — z w z + Azw + — (m + n + ran — m — n) = 0. (6-1) 
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Figure 3: Configuration of zeros for the Ohyama polynomial Pj(z) 

The Laurent expansion for solutions of this equation near movable singularity 
z — zq is 

f \ 1 Zq 2 / \ ^0 / \2 

w(z) = \-ai — (Z - Zq) + — [Z- Zq) + 

Z — Zq 9 18 



2z ai 4 4 1 4 4 4 2 

— ran H m zn H n 

' - 45 45 



+ — n m 

45 45 15 45 45 405 



(6.2) 



(z - Zq) 3 + ( y - — ) :n i ' «<; U 



5 , Zq 



9720 



[8 2o 4 - 



-432 z ai + 288 mn - 288 n + 288 m 2 - 363 + 288 n 2 - 288 to] (z - z ) 

Then relations for the zeros of the generalized Okamoto polynomials can be 
written as 



S (zfc ~^ )+fll °' 



(6.3) 



E 



^ (z k - z 3 ) 2 9 



-4 = o, 



(6.4) 



12 



E 



I + 51 = 



(6.5) 



E 



405 



2 ai 2j 
15 



(6.6) 



H (l + m + rt — m 2 — nin — n 2 ) =0, 

45 v ; 

Okamoto polynomial P n ,m{z) 



m A + n 2 + mn 
The generalized 



of the generalized 



where the degree p 

1 i 2 i 

. .- i ... i "™, — TTj, — n. 

. Okamoto polynomials Pa^z) takes the form [12] 



is p 



.40 



1755 o fi 4238325 310134825 2S 
— — z 36 H — z 32 z 28 - 



359947121625 24 40989148261425 20 5340997262498625 16 

128 Z 256 Z + 2048 * 

299607997122354375 12 35054135663315461875 8 

2048 " 65536 

849388671841874653125 4 1528899609315374375625 

131072 Z + 1048576 



+ - 



(6.7) 
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Figure 4: Configuration of zeros for the Okamoto polynomial P 4j 4(z) 

The configuration of zeros for the Okamoto polynomial P 4j 4(z) are demon- 
strated in the Figure 4. 
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7 Relations for the zeros of the generalized Her- 
mite polynomials 

Applying the transformation y z — w{z)y we get from (|1.6|l the following third- 
order differential equation 

w zzz + 6 w z 2 — 4 w z z 2 — 8 w z (n — m) + 4 z w — 8 m n = 0. (7-1) 
The power expansion for its solutions near the point z = zq is 

uj(z) = - — - — - + ai + \ (z 2 + 2 n - 2 to) (z - z ) + - z (z - z ) 2 + 

+ -^r ^3 - nzo 2 - n 2 - TOn + toz 2 - to 2 + ^z ai - ^ Z Q 4 ^ (-2 - zof + (7-2) 



(2toz - z 3 - 2?iz + 3 ax) (z - z ) i + a 6 (z - z )° + ... 

18 

Using the series l|7.2[) we obtain the following relations for the zeros of the 
generalized Hcrmitc polynomials 



nm 



(7.3) 



nri i 



v- 1 z 2 2(n-m) 



nm ^ 

v ~ ^ - 0, (7.5) 



(*-*>)» 2 



nri i 



y-v 1 zj ^4(m — n)zj ^ 2 ai ^ 



(7.6) 



^^.(zfc-z,) 4 45 45 15 



H (3 — to 2 — inn — n 2 ) = 0. 

45 v 1 

The generalized Hermite polynomial H^(z) is [14] 



lM,)_,» + ^«-f* + ™^ + »» (7.7, 

The configuration of the zeros of the generalized Hermite polynomial H^^z) 
are illustrated in the Figure 5. 
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Figure 5: Configuration of zeros of the generalized Hermite polynomial H^(z) 



8 Conclusion 

In this work we have studied special polynomials associated with the rational 
solutions of the second, the third, and the fourth Painleve equations. We have 
shown that if a family of polynomials satisfies an n-linear differential equation, 
then much information concerning the polynomials can be obtained from the 
asymptotic study of the solutions of a differential equation related to the original 
n-lincar one. In particular, it can be found the degree of each polynomial, 
symmetric functions of its roots and some other correlations for the roots called 
the generalized Stieltjes relations. Also one can obtain possible degrees of the 
roots. In this work we have derived relations for the roots of the following 
polynomials: the Laguerre polynomials, the Umemura polynomials, the Ohyama 
polynomials, the generalized Okamoto polynomials, and the generalized Hermite 
polynomials. Our approach can be applied to other families of polynomials if 
there exists a differential equation satisfied by the family in question [28-31]. 
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